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^  ABSTRACT 

A  quantitative  model  Is  developed  for  computing  currents  and  electric 
fields  In  the  area  covered  by  the  poleward  set  of  field-aligned  currents 
(region-1  Birkeland  currents).  The  high-latitude  band  considered  carries  most 
of  the  westward  electrojet  and  lies  just  poleward  of  the  region  treated  in  the 
Rice  University  computer  simulation  (which  covers  the  inner  magnetosphere  and 
corresponding  ionosphere);  the  present  high-latitude-current  model  supplements 
the  inner-magnetosphere  simulations  to  make  them  nearly  global.  Birkeland 
current  and  auroral ly  enhanced  conductivity  are  assumed  to  be  uniformly  dis¬ 
tributed  across  the  band.  The  area  poleward  of  the  band  is  taken  to  be  an 
Insulator,  while  boundary  conditions  at  the  low-latitude  edge  of  the  band  are 
derived  from  the  results  of  the  lower-latitude  computer  simulations  of  the 
substorm  event  that  occurred  on  September  19,  1976.  The  time-dependent  con¬ 
ductivity  model  used  is  based  on  electron  fluxes  and  mean  energies  measured 
from  the  Air  Force  S3-2  satellite.  Ionospheric  electric  fields  and  currents 
are  calculated;  ionospheric  neutral  winds  are  neglected. 

Model  predictions  of  the  locations  of  the  electric-field  reversal  agree 
well  with  typicaul  0G0-6  satellite  observations^  .The  model  also  agrees  with 
S3-2  data  from  September  19,  1976  on  an  important  feature,  namely,  the  elec¬ 
tric  field  reversal  usually  occurs  poleward  of  the  bulk  of  region-1  current. 

y,  Joule  heating  of  the  upper  atmosphere,  as  inferred  from  calculation  of 
currents  and  electric  fields  in  the  high-latitude  region  considered  here,  was 
found  to  be  about  2  x  watts  during  the  substorm  period,  somewhat  less 
than  the  value  estimated  for  the  lower-latitude  region  (region-2  and  below). 
Model  values  for  the  strength  of  the  electrojet  and  the  amount  of  Joule  heat¬ 
ing  agree  to  within  ~  20X  with  results  based  on  a  simple  Cowling  conductivity 
band  (with  no  latitudinal  current  flowing  into  the  band). 
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I.  INTRODUCTION 


Work  at  Rice  University  over  the  past  few  years  has  produced  a  method  of 
computer  simulating  convection  in  the  inner  magnetosphere  (L  <  10).  The  Rice 
Convection  Model  self-consistently  includes  ionospheric,  Birkeland  and  ring 
currents  [Jaggi  and  Wolf,  1973;  Harel  and  Wolf,  1976;  Harel  et  al.,  1981a). 
Recently,  two  geophysical  events — the  substorm-type  event  of  September  19, 
1976  and  the  magnetic  storm  of  July  29,  1977 — have  been  computer  simulated, 
and  theoretical  predictions  have  been  compared  with  a  variety  of  observational 
data  [Harel  et  al.,  1981b;  Spiro  et  al.,  1981;  Wolf  et  al.,  1981).  The  region 
of  applicability  of  the  Rice  Convection  Model  has  been  restricted  to  the  inner 
magnetosphere.  Uncertainties  about  the  detailed  physics  governing  solar-wind/ 
magnetosphere  coupling  and  magnetotail  dynamics  have  prevented  us  from 
including  the  magnetopause  boundary  layer,  the  tail  lobe  and  the  outer  plasma 
sheet  in  the  self-consistent  modeling  scheme.  In  other  words,  we  exclude  from 
the  self-consistent  simulations  [Harel  et  al.,  1981a)  the  part  of  the  magneto¬ 
sphere  that  serves  as  the  generator  for  magnetospheric  convection  [Stern, 
1977).  Part  of  the  auroral  ionosphere  has  also  been  excluded,  namely,  the 
polar  cap  and  the  part  of  the  ionosphere  that  corresponds  to  the  region- 1 
Birkeland  currents,  which  connect  the  ionosphere  to  the  generator  of  the 
magnetospheric-convection  system.  [See,  e.g.,  Schield  et  al.,  1969;  Wolf, 
1974;  or  Stern,  1977. J  The  effect  of  the  generating  currents  on  inner-magne- 
tospheric  convection  has  been  represented  by  a  boundary  condition. 

Exclusion  of  the  high-latltude  part  of  the  auroral  zone  from  the  modeling 
region  has  proved  awkward  in  several  respects.  Most  of  the  westward  electro¬ 
jet  flows  in  the  high-latitude  part  of  the  auroral  zone.  Neglecting  the  high- 
latitude  ionospheric  currents  consequently  precludes  any  serious  comparisons 
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between  theoretical  and  observed  ground-magnetometer  data.  It  also  precludes 
accurate  calculation  of  global  Joule  heating. 

The  present  model  represents  an  extension  of  the  Rice  Convection  Model  to 
the  highest-latltude  part  of  the  auroral  ionosphere.  Although  we  cannot  yet 
reliably  computer  simulate  the  physics  of  the  far-magnetosphere  ends  of  these 
field  lines,  we  can  model  their  ionospheric  ends,  using  observational  infor¬ 
mation  and  simple  assumptions  about  the  distribution  of  current  coming  down 
from  the  magnetosphere. 

Many  other  theoretical  models  [Yasuhara  et  al.,  1975;  Yasuhara  and 
Akasofu,  1977;  Gizler  et  al.,  1979;  Nisbet  et  al.,  1978;  Nopper  and  Caro- 
villano,  1978,  1979;  Kamide  and  Matsushita,  1979a,  1979b J  have  treated  the 
relationship  between  horizontal  ionospheric  currents  and  field-aligned  cur¬ 
rents.  However,  in  these  approaches,  "typical"  or  "average”  conditions  of 
magnetic  activity  were  considered,  and  only  a  general  comparison  with  data  was 
possible.  Most  studies  use  an  average  Birkeland  current  distribution  as 
input,  with  a  relatively  simple  conductivity  model.  Akasofu  et  al.  [1980] 
have  used  data  from  the  Alaska  meridian  magnetometer  chain  to  infer  Birkeland 
and  horizontal  ionospheric  currents.  On  the  other  hand,  the  Rice  Convection 
Model  [Harel  et  al . ,  1981a,  1981b]  does  detailed  simulations  of  individual 
events  for  more  precise  comparison  with  various  types  of  data. 

The  model  described  In  this  paper  complements  the  Rice  Convection  Model. 
The  model  uses  polar  cap  potential  drop  and  fairly  realistic  conductivity  as 
input  to  compute  Birkeland  currents,  electric  fields,  and  horizontal  iono¬ 
spheric  currents. 

The  modeled  region  is  a  band  ~  3°  wide,  centered  at  about  70-73°  in  geo¬ 
magnetic  latitude.  More  precisely,  the  low-latitude  boundary  of  the  band  dis¬ 
cussed  in  this  paper  is  the  low-latitude  edge  of  region-l  Birkeland  current. 
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The  poleward  boundary  of  Che  band  is  defined  by  the  poleward  extent  of  the 
aurorally  enhanced  conductivity,  as  inferred  from  satellite  observations.  The 
exact  locations  of  the  boundaries  are  variable. 

The  following  section  presents  a  detailed  explanation  of  the  theory.  In 
the  third  section,  this  theory  is  applied  to  a  specific  substorra-type  event 
that  occurred  on  September  19,  1976.  East-west  "electrojet"  currents  and 
Joule  heating  rates  are  calculated  for  our  modeling  region  for  this  event. 
Birkeland  current  strengths  and  electric  field  reversal  locations  are  compared 
with  average  or  typical  satellite  data.  The  impact  of  these  high-latitude 
currents  on  ground  magnetic  variations  is  discussed  in  a  separate  paper  [Chen 
et  al.,  1981]. 

II.  THEORY 


Current  conservation  implies  that  field-aligned  current  is  balanced  by 
the  divergence  of  horizontal  ionospheric  current.  Thus,  current  conservation 
can  be  expressed  as 


is  used  (neglecting  neutral  winds),  where  E  is  the  electric  field  in  the  frame 
of  reference  that  rotates  with  the  solid  earth,  and  E  is  the  tensor  repre¬ 
senting  height-integrated  conductivity  [Wolf,  1970;  Jaggi  and  Wolf,  1973; 
Rishbeth  and  Garriott,  1969,  p.  138]: 

E-(^  »> 

*5  ** 

Here,  ZCe  =  Ep/sin2  I,  =  -E^  *  E^/sin  I,  and  E^  =  Ep,  where  Ep  and  E^ 

are  height-integrated  Pedersen  and  Hall  conductivities,  respectively,  and  I  is 
the  magnetic  dip  angle.  Combining  (1),  (2),  and  E  =  -  7  V,  where  V  is  the 
electrostatic  potential  (assuming  that  any  induction  fields  are  negligible), 
we  get: 


7  •  (J  •  V  V)  =  j„.  (4) 

Within  the  band,  conductivity  depends  on  local  time,  but  is  independent  of 
latitude.  A  thin  band  is  assumed  (i.e.,  the  band  width  is  small  compared  to 
the  distance  to  the  pole).  For  a  thin  band,  the  azimuthal  component  of 
electric  field  is  approximately  independent  of  distance  across  the  band, 
since  7  x  E  =  0  implies  (by  Stokes'  law)  that  J  E  •  dl  m  0.  We  assume 
Birkeland  current  to  be  independent  of  latitude  within  the  modeling  region. 
With  this  assumption,  and  within  the  thin  band  approximation,  current  conser¬ 
vation  implies  that  horizontal  ionospheric  currents  and  the  poleward  component 
of  the  electric  field  vary  almost  linearly  with  latitude. 

Figure  l  presents  a  definition  of  coordinates  used  in  the  model.  As  an 
approximation,  the  polar  cap  is  taken  to  be  flat.  A  two-dimensional  polar 
coordinate  system  is  used.  (The  center  of  the  coordinate  system  is  offset 
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from  the  geomagnetic  pole  by  two  degrees  toward  midnight).  Boundary  condi¬ 
tions  are: 


V  =  V (\p )  at  f  «  a, 

(5) 

jf  =  SW  at  £  *  a,  and 

(6) 

=  0  at  5  =  b(4>). 

(7) 

(Equation  (7)  follows  from  the  approximation  that  the  polar  cap  is  an  insu¬ 
lator.)  We  write  the  azimuthal  electric  field  in  the  form: 


P  (r  -  -  1  3V(C.*)_  1  dV(*) 


(8) 


(In  the  thin-band  approximation,  we  neglect  the  independence  of  3V(£  ,i|>)/3i|>. ) 
The  latitudinal  components  of  electric  field  at  b  and  a  are  evaluated  using 
equation  8,  the  boundary  conditions  given  by  equations  6  and  7,  and  Ohm's  law 
for  ionospheric  currents  (equation  2).  We  obtain: 


EC<M) 


E^(a,t|») 


lg(<V)  +  £ 


(♦>  i 

a  dip  1 


(9) 


(10) 


We  assume  that  E^(^,t|/)  varies  linearly  between  b  and  a.  Using  the  linear 
dependence  of  Ej-  on  %  ,  and  integrating  leads  immediately  to  an  expression  for 
V(C  ,<!<)•  Evaluation  of  this  expression  at  £  *  a  confirms  self-consistency  of 
the  potential  used  with  assumed  boundary  conditions. 

Prom  a  knowledge  of  E  (£  ,\|>)  and  E (£,i|/),  total  east-west  current  (i.e., 

T 

the  "electrojet"),  and  Joule  heating  can  be  calculated.  Total  eastward  cur¬ 
rent  is  given  by: 
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\  -  /bV£  •  V*’  -i*1  ‘"(j>  -  V*)(“ - 5-1 - )<•  - b>-  <“> 


Joule  heating  (per  unit  ip)  is  given  by: 


JH  =  /Jj  •  =  /S{EK(*)[Ee(5,*)]S 


+  E^<*)l%<e.<»)]2}cdc. 


This  becomes : 


JH  -  tM  «,)  (^<£1)2  in  <./b)  +  £?{  <*)  (^)([E?(b,*>]2  (i+Jt) 


12  r  a  +  3b>| 


+  Ef(b,*)  E  (a,*)(a  +  b)  +  [E  (a,*)]*  (t+Ji)}. 


With  ionospheric  current,  Birkeland  current  is  calculated  from  current 
conservation.  Multiplying  (1)  by  Sd£  and  integrating  from  £  *  b  to  £  *  a 
gives 


r»  <  ,*  .  r«  8({V 

'  b  h£d£  Jb  3E 


This  is  equivalent  to: 


dJ  dJ,  dJ_ 

- 1  a - 1.  JL 

dip  dip  dip  ’ 


dJ,  dJr 

Ii  r 

where  and  — -  are,  respectively,  the  total  Birkeland  and  equatorward  cur¬ 
rents  leaving  the  band  per  unit  ip. 
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III.  APPLICATION  OF  MODEL  TO  SEPTEMBER  19,  1976  SUBSTORM  EVENT 


Harel  et  al.  [1981a,  1981b]  have  simulated  the  substorm  event  of  Septem¬ 
ber  19,  1976.  Figure  3  of  Harel  et  al.  [1981a]  shows  the  Fort  Churchill  H 
magnetogram  and  the  estimated  potential  drop  for  the  event.  Onset  of  the 
event  occurred  at  approximately  1000  UT.  Figures  2a  and  2b  show  the  equator- 
ward  boundary  conditions  that  were  furnished  by  the  original  simulation  and 
were  used  as  input  to  the  present  hlgh-latitude  model.  The  figures  correspond 
to  1100  UT,  near  the  observed  peak  of  substorra  activity. 

The  conductivity  model  considers  auroral  enhancement  and  photoionization. 
(Our  method  of  including  the  sunlight  effect  Is  described  by  Harel  et  al. 
[1981a,  pg.  2234 ] - )  At  the  high  latitudes  considered  in  the  present  model, 
auroral  enhancement  is  the  dominant  effect.  We  use  equations  (9)  and  (10)  of 
Harel  et  al.  [1981a]  to  infer  auroral  conductivity  enhancements  from  electron 
fluxes  and  average  energies  measured  by  the  S3-2  satellite  [D.  A.  Hardy,  pri¬ 
vate  communication,  1977],  (See  Wallis  and  Budzlnski  [1981]  for  a  comparison 
of  equations  [9]  and  [10]  with  more  precise  calculations.)  Integration  of 
satellite  data  over  a  distance  defined  by  electron  flux  enhancement,  and  then 
fitting  the  data  to  a  trigonometric  function  of  local  time,  provides  the 
auroral  enhancement  conductivity  model.  Actual  measurements  were  near  dawn 
and  dusk,  so  the  that  all  other  local  times  were  obtained  by  extrapolation. 
The  model  allows  auroral  conductivity  enhancement  at  midnight  to  be  ~  4  times 
greater  than  the  conductivity  enhancement  at  noon,  for  the  high-latitude  area 
treated  here.  To  obtain  the  total  conductivity,  contributions  from  the 
auroral  enhancement,  photoionization,  and  base-level  nightside  values  are 
summed  according  to  equations  (A. 20)  and  (A. 21)  of  Harel  et  al.  [1981a]. 
Figure  2c  depicts  the  local-time  dependent  conductivity  for  1100  UT.  The 
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width  of  the  band  also  depends  on  electron-flux  measurements.  Integration  of 
satellite  data  over  a  distance  defined  by  the  region  of  high  electron  flux, 
and  then  fitting  the  data  to  a  trigonometric  function  of  local  time,  provides 
the  model  for  the  width  of  the  band.  The  poleward  boundary  is  taken  to  be  a 
function  of  local  time,  but  the  equatorward  boundary  remains  constant  in  this 
approach. 

The  Birkeland  current  pattern  predicted  from  this  model  is  compared  with 
the  observational  data  of  Iijima  and  Potemra  [1979,  private  communication)  in 
Figure  3.  Included  in  this  figure  is  a  plot  (from  Figure  2b)  of  -j^(a)  (which 
is  equal  to  Birkeland  current  if  one  neglects  east-west  current  for  £  <  a). 
The  present  calculation  improves  upon  this  approximation  by  including  effects 
of  the  high-latitude  auroral  electro jet  (which  is  known  to  intensify  during  a 
substorm).  It  should  be  emphasized  that  plotted  data  were  statistically  aver¬ 
aged  for  |AL|  >  100  y  by  Iijima  and  Potemra.  The  model  predicted  Birkeland- 
current  strengths  tend  to  be  larger  than  the  observational  average,  probably 
because  |AL|  was  about  600  y  at  1100  UT. 

An  analysis  of  the  dependences  of  Birkeland  current  and  the  electrojet 
(within  the  modeled  region)  on  boundary  conditions  provides  important  physical 
Insights.  In  this  model,  Birkeland  current  is  mainly  “controlled"  by  the 
latitudinal  current  at  the  equatorward  boundary.  Additionally,  in  a  simple 
Cowling  model,  with  j^  =  0  everywhere  within  the  thin  band,  we  have 


so  that  the  electrojet  is  equal  to  the  product  of  the  Cowling  conductivity, 
the  azimuthal  component  of  electric  field  and  the  band  width.  Figure  4  shows 
that  values  from  the  Cowling  conductivity  picture  agree  with  the  present  model 
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to  within  ~  20%.  Apparently,  the  total  east-west  current  is  mainly  "con¬ 
trolled"  by  east-west  electric  field,  which  was  essentially  a  boundary  condi¬ 
tion. 

As  illustrated  in  Figure  4,  the  calculated  westward  electrojet  in  the 
region  treated  in  this  paper  (from  equation  11)  extends  from  about  1900  LT 
past  midnight  to  about  0900  LT,  with  a  maximum  strength  of  ~  11  x  105  amp  in 
post-midnight  hours.  (This  is  in  reasonable  agreement  with  observational  data 
for  substorm  conditions.)  Figure  5  shows  total  ionospheric  current  density 
for  geomagnetic  latitudes  poleward  of  ~  60°.  Current  vectors  displayed 
include  both  results  from  the  present  model  and  the  lower  latitude  Rice  simu¬ 
lation.  Most  of  the  westward  electrojet  lies  in  the  high-latitude  area  of 
region-1  current.  The  calculated  eastward  electrojet  is  comparable  to  (but 
slightly  weaker  than)  the  portion  of  the  eastward  electrojet  in  the  lower- 
latitude  region;  the  lower  latitude  electrojet  extends  to  ~  2300  LT. 

Theoretical  magnetograms  were  computed  and  compared  with  observational 
data  for  the  substorm  event  of  September  19,  1976.  Magnetic  perturbations 
observed  on  the  ground  result  mostly  from  the  current  systems  in  the  Earth's 
ionosphere  and  magnetosphere.  The  details  of  the  method  used  to  calculate 
theoretical  magnetograms  can  be  found  in  Chen  et  al.  [1981].  The  Birkeland 
current  and  east-west  currents  in  region-1  are  part  of  a  larger  current 
system.  It  should  be  emphasized  that  since  the  westward  electrojet  flows 
mainly  in  the  poleward  region  and  has  large  ground  magnetic  effects,  the  addi¬ 
tion  of  the  auroral  band  described  in  this  paper  greatly  improves  agreement 
between  theoretical  and  observed  magnetograms.  The  comparison  of  theoretical 
and  observed  magnetograms  is  presented  and  discussed  by  Chen  et  al.  [1981). 

In  the  present  model,  it  also  is  possible  to  Investigate  the  location  of 
electric  field  reversals.  Figure  6  presents  a  diagram  of  the  thin  band  and 
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the  predicted  global  electric  field  reversal  pattern  for  geomagnetic  latitudes 
poleward  of  ~  50°.  The  global  pattern  predicted  from  the  model  is  in  very 
good  agreement  with  the  observed  pattern  shown  in  Figure  7.  In  both  Figures  6 
and  7,  the  reversal  is  at  ~  75°  geomagnetic  latitude  near  dusk  and  dawn, 
whereas  near  midnight  (in  the  classical  Harang  discontinuity  region),  the 
reversal  appears  at  ~  63°  geomagnetic  latitude.  Overall  concurrence  of  Figure 
6  with  observations  is  especially  encouraging. 

Values  for  Joule  heating  derived  from  this  model  (from  equation  13)  are 
shown  graphically  in  Figure  8a  for  UT  =  1100.  In  the  simple  Cowling  conducti¬ 
vity  example  discussed  previously.  Joule  heating  per  area  is  given  by  the 
expression: 


JH 


(16) 


which  is  equivalent  to 


JH  = 


(17) 


In  this  picture.  Joule  heating  (per  length)  is  equal  to  the  product  of  the 
Cowling  conductivity,  the  square  of  the  azimuthal  component  of  the  electric 
field,  and  the  band  width.  Figure  8a  shows  that  values  from  the  Cowling  con¬ 
ductivity  picture  (with  no  latitudinal  current  flowing  into  the  band)  agree 
with  the  present  model  to  within  ~  20%.  Since  the  azimuthal  component  of 
electric  field  is  proportional  to  the  derivative  of  the  equatorward  boundary 
potential  (equation  8),  the  above  discussion  implies  that  Joule  heating  in  the 
high-lat itude  band  is  mainly  "controlled"  by  the  equatorward  boundary  poten¬ 
tial.  A  detailed  background  discussion  on  Joule  heating  is  provided  in  Harel 
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et  al.  [1981b],  along  with  results  of  their  study.  This  effort  helps  to  make 
the  previous  investigation  more  global.  Figure  8b  shows  Joule  heating  contri¬ 
butions  from  both  the  present  high-latitude  model  (indicated  by  circles)  and 
the  lower  latitude  simulation  for  1100  UT.  The  high-  and  low-latitude  regions 
tend  to  complement  each  other,  so  that  the  total  Joule  heating  is  relataively 
uniform  as  a  function  of  local  time.  Included  in  Figure  9  are  graphs  indi¬ 
cating  Joule  heating  and  ring  current  energy  as  functions  of  universal  time. 
The  top  panel  indicates  that  the  integrated  (total)  Joule  heating  is  about  3.2 
times  the  change  in  ring  current  energy  for  the  modeled  substorm.  (See  Harel 
et  al.  [1981b].)  The  bottom  panel  shows  that  Joule  heating  poleward  of  the 
equatorward  edge  of  region-1  currents  comprises  ~  40%  of  the  total  Joule  heat¬ 
ing  in  our  modeled  event.  Clearly,  the  region  we  treat  here  makes  a  signifi¬ 
cant  contribution  to  global  Joule  heating. 

The  prime  incentive  for  this  study  was  to  investigate  the  extent  to  which 
the  higher  latitude  area  (poleward  of  the  main  simulation  region  of  the  Rice 
Convection  Model)  affected  ground  magnetic  variations  and  Joule  heating  of  the 
upper  atmosphere.  Addition  of  this  auroral  band  of  high  conductivity  has 
significantly  improved  calculated  magnetograras.  Improved  agreement  between 
observed  and  calculated  magnetograms  is  due  to  inclusion  of  the  large,  often 
dominant,  ground  magnetic  effects  of  the  westward  electrojet  flowing  mainly  in 
the  higher  latitude  region  treated  in  this  paper.  In  addition,  this  investi¬ 
gation  Indicated  that  the  amount  of  Joule  heating  of  the  upper  atmosphere 
inferred  from  calculations  of  currents  and  electric  fields  in  the  higher  lati¬ 
tude  region  is  significant. 
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FIGURE  CAPTIONS 


Figure  1. 


Figure  2. 


Figure  3. 


Figure  4. 


Polar  view,  illustrating  definitions  of  coordinates  used  in  the 
model . 

(a)  Equatorward  boundary  potential  (used  as  input  to  model)  as 
a  function  of  local  time  (UT  *  1100). 

(b)  Co-latitudinal  component  of  horizontal  ionospheric  current 
at  the  equatorward  boundary  (used  as  input  to  model)  as  a  func¬ 
tion  of  local  time  (UT  ■  1100). 

(c)  Conductivity  (used  as  input  to  model)  as  a  function  of 
local  time  (UT  *  1100).  ["o'*  refers  to  ,  ”x"  refers 
to  ”  "  refers  to  S^.J 

Birkeland  current  as  a  function  of  local  time.  Calculations 
from  this  paper  are  denoted  by  (UT  =  1100).  For  the  simpler 
picture,  in  which  *  -j^  (a),  the  curve  is  denoted  by  "x”. 
Data  from  Iijima  and  Potemra  (1979,  private  communication)  are 
indicated  by  "A”. 

Total  east-west  current  in  our  modeling  region  as  a  function  of 
local  time.  Calculations  using  boundary  conditions  from  lower 
latitude  modeling  are  denoted  by  "o".  Calculations  using 
V(a)  =  0  are  denoted  by  "x".  Calculations  using  j^(a)  -  0 
(Cowling  conductivity  picture)  are  denoted  by  "A"  (UT  «  1100). 
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Figure  5 


Figure  6. 


Figure  7. 


Ionospheric  total  current  density.  The  two  highest  latitude 
"circles"  of  arrows  are  at  the  poleward  and  equatorward  bound¬ 
aries  of  the  present  high-latitude  model.  Arrows  at  lower  lati¬ 
tudes  are  from  results  of  the  lower  latitude  computer  simulation 
(UT  »  1050). 

Diagram  of  places  where  polarity  reversals  of  radial  electric 
field  occur  for  UT  «  1100.  The  dotted  region  is  the  band 
modeled  here.  The  cross  represents  the  geomagnetic  north  pole. 
The  small  dot  at  the  center  is  the  2°  offset  discussed  in  the 
text.  refers  to  a  reversal  being  inside  the  band.  "x" 
refers  either  to  a  predicted  reversal  being  poleward  of  the 
region  studied  (which  physically  means  that  the  reversal  must 
occur  on  the  poleward  boundary,  as  shown),  or  to  a  predicted 
reversal  being  at  the  equatorward  boundary.  "A"  refers  to 
predicted  reversal  positions  from  lower  latitude  (to  50°) 
studies  of  Harel  et  al.  (1981a, b). 

Schematic  representations  of  typical  convective  flow  and  elec¬ 
tric  field  directions  looking  down  on  north  magnetic  pole.  The 
boundaries  are  for  moderately  disturbed  to  disturbed  conditions. 
From  Fig.  1  of  Maynard  [1974]. 
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Figure  8. 


(a)  High- latitude  Joule  heating  per  unit  (both  hemispheres) 

as  a  function  of  local  time.  (See  equations  12  and  13.)  "o" 

refers  to  calculations  using  boundary  conditions  from  lower 
latitude  modeling;  "x“  refers  to  calculations  for  V(a)  =  0;  "A” 
refers  to  calculations  for  j^(a)  =  0,  i.e.,  the  Cowling-conduc¬ 
tivity  picture  (UT  =  1100). 

(b)  Joule  heating  per  unit  (both  hemispheres)  as  a  function 
of  local  time  showing  contributions  from  both  the  high  and  low 
latitude  models.  "o"  designates  calculations  from  the  present 
high-latitude  model.  Results  from  the  lower  latitude  computer 
simulation  are  Indicated  by  a  solid  line,  with  no  symbols  or 
label  attached.  The  topmost  curve  is  the  sum  of  calculations 
from  both  the  high  latitude  model  and  lower  latitude  simulation 
(UT  -  1100). 
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Figure  9 


The  top  panel  shows  ring  current  energy  and  integrated  (total) 
Joule  heating  as  functions  of  universal  time  through  the  event. 
The  bottom  panel  shows  Joule  heating  (both  hemispheres)  as  a 
function  of  universal  time  through  the  modeled  event.  Calcula¬ 
tions  from  this  paper,  i.e.,  due  to  the  region  of  the  ionosphere 
poleward  of  the  equatorward  edge  of  the  region-1  currents,  are 
indicated  by  "x".  Calculations  from  Harel  et  al.  [1981b], 
covering  the  region  of  the  ionosphere  equatorward  of  the  equa¬ 
torward  edge  of  region-1  currents,  are  indicated  by  a  solid 
line,  with  no  symbols  or  label  attached.  The  topmost  curve  is 
the  sum  of  the  other  two  curves,  i.e.,  total  Joule  heating  cal¬ 
culated  in  both  the  present  high-latitude-current  model  and  the 
lower  latitude  computer  simulation.  Initial  substorm  onset  was 
at  10  UT,  and  activity  remained  high  through  12  UT. 
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